C||  A(l'Oo<}-66o 

£. 


F-1 1 1C  Longitudinal  and  Lateral 
Aerodynamic  Flight  Data  Analysis  for 
Take-off  and  Landing  Configurations 

A.D.  Snowden  and  J.S.  Drobik 


DEFFNCH  SCIENCE  ANi)  TEClINCiJXn  ORGAMSATiON 


F-lllC  Longitudinal  and  Lateral  Aerodynamic  Flight  Data 
Analysis  for  Take-Off  And  Landing  Configurations 


DSTO-TR-0321 

Andrew  D.  Snowden  and  Jan  S.  Drobik 

Air  Operations  Division 
Aeronautical  and  Maritime  Research  Laboratory 


ABSTRACT 

A  series  of  flight  trials  was  performed  on  an  F-lllC  aircraft,  tail  number  A8-132,  at  the 
RAAF's  Aircraft  Research  and  Development  Unit.  Data  obtained  from  the  tests  has  been  anal¬ 
ysed  in  Air  Operations  Division  to  determine  the  aircraft  aerodynamic  and  control  derivatives. 
In  this  report  the  longitudinal  and  lateral  derivatives  for  the  take-off  and  landing  configurations 
are  presented  and  compared  with  wind  tunnel  results. 
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F-lllC  Longitudinal  and  Lateral  Aerodynamic  Flight  Data 
Analysis  for  Take-Off  And  Landing  Configurations 

EXECUTIVE  SUMMARY 

With  the  introduction  of  the  variable  wing  sweep  General  Dynamics  F-lllC  into  service  in  1973, 
the  Royal  Australian  Air  Force  (RAAF)  was  faced  with  providing  technical  support  for  a  range 
of  complex  systems.  In  particular,  the  aircraft  has  a  complex  adaptive  flight  control  system  which 
is  needed  to  provide  satisfactory  flying  behaviour  across  the  aircraft's  large  operational  flight 
envelope.  At  the  time  of  acquisition  of  the  F-lllC,  it  was  not  possible  to  obtain  from  the  manu¬ 
facturers  a  flight  dynamic  model  of  the  aircraft  suitable  for  flight  d3mamic  analysis.  To  provide 
these  analysis  capabilities  a  model  was  developed  at  the  Air  Operations  Division  (AOD)  based 
on  the  manufacturers  design  data  reports  and  from  wind  tunnel  tests  in  the  AOD  Transonic 
Wind  Tunnel.  The  model  enables  the  investigation  of  the  effects  of  pilot  commands,  external 
stores,  control  system  modifications  and  atmospheric  turbulence  on  the  rigid  body  response. 

The  F-lllC  aircraft  operated  by  the  RAAF  has  a  different  geometry  from  the  versions  oper¬ 
ated  by  the  United  States  Air  Force  (USAF).  The  aerodynamic  data  supplied  to  Australia  with 
the  F-lllC  included  data  developed  for  the  F-lllA  aircraft  which  has  the  same  fuselage,  and  for 
the  F-lllB  which  has  the  same  wing.  To  validate  these  aerodynamic  data,  a  flight  test  program 
was  conducted  at  the  RAAF's  Aircraft  Research  and  Development  Unit  (ARDU)  to  record  the 
flight  behaviour  of  an  F-lllC  in  manoeuvres  designed  to  optimise  the  identifiability  of  aero¬ 
dynamic  characteristics.  Analysis  of  the  data  has  been  carried  out  by  AOD  to  determine  the 
aerodynamic  characteristics. 

Data  from  the  program  has  been  used  to  validate  a  comprehensive  flight  dynamic  model  of 
the  F-lllC  which  has  been  developed  by  AOD.  Data  from  the  model  is  currently  being  used  in 
support  of  the  F-lllC  Mission  Simulator  and  Avionics  Update  programs  to  enhance  the  opera¬ 
tional  effectiveness  of  the  aircraft. 

Stability  and  control  derivatives  describing  the  longitudinal  and  lateral  aerodynamic  char¬ 
acteristics  of  the  F-lllC  have  been  determined  from  flight  test  measurements.  The  results  pre¬ 
sented  in  this  report  are  for  16°  and  26°  wing  sweep  angles  in  the  take-off  and  landing  configu¬ 
ration. 

In  general,  all  primary  longitudinal  and  lateral  derivatives  show  satisfactory  repeatability 
and  trends  like  those  of  the  General  Dynamics  and  AOD  wind  tunnel  data.  However,  the  air¬ 
craft  pitch  stiffness,  yaw  stiffness  and  rudder  effectiveness  were  consistently  lower  in  magni¬ 
tude  than  the  corresponding  wind  tunnel  data  for  both  wing  sweeps.  This  is  consistent  with 
previous  results  for  the  clean  aircraft  configuration. 
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NOTATION 


Normal,  longitudinal  and  lateral  acceleration  (g) 
Angle-of-attack  and  sideslip  measurement  bias  (°) 

Axial  force  coefficient 
Rolling  moment  coefficient 
Pitching  moment  coefficient 
Normal  force  coefficient 
Yawing  moment  coefficient 
Longitudinal  force  coefficient 
Side  force  coefficient 
Reference  chord  (ft) 

Centre  of  gravity 
Gravitational  acceleration 
Altitude  (ft) 

Moments  of  inertia  about  roll,  pitch  and  yaw  axes  (slug.ft^) 
Cross  product  of  moment  of  inertia  (slug.ft^) 

Interval  counter 
Cost  function 

Flow  amplification  factors  for  angle-of-attack  and  sideslip 
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Degrees  per  radian  (57.2958) 

Yaw  rate  (°/s) 

Reference  area  (ft^) 
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NOTATION  (continued) 
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A 
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Software  Variables 

Quantities  listed  in  parenthesis  indicate  software  notation  used  in  time  history  figures  in  this 
report. 

Subscripts 

CADS  Central  air  data  system  instrument 

CM  Crew  module  instrument 

eg  Centre  of  gravity 

i  Interval  counter 

L  Left  (port) 

m  Measured  quantity 

NBTU  Nose  boom  transducing  unit  instrument 

p,  q,  r  Rate  derivatives  with  respect  to  indicated  quantity 

R  Right  (starboard) 

true  True  value 

z  Computed  response 
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Thrust  (Ibf) 

Time  (s) 

Control  vector 
Airspeed  (kn) 

Response  weighting  matrix 
State  vector 

Longitudinal  instrument  offsets  from  eg  (ft) 

Lateral  instrument  offsets  from  eg  (ft) 

Measured  response  vector 
Computed  response  vector 
Vertical  instrument  offsets  from  eg  (ft) 

(ALPHA)  Angle-of-attack  (°) 

(BETA)  Angle-of-sideslip  (°) 

Control  deflection  (°) 

Aileron  deflection  i.e.  differential  stabilator  -  (°) 

Elevator  deflection  i.e.  symmetrical  stabilator  (°) 

Flap  deflection  (°) 

Rudder  deflection  (°) 

Spoiler  deflection  (°) 

Stabilator  deflection  (°) 

(THETA)  Pitch  angle  (°) 

Wing  sweep  angle  (°) 

Angle-of-attack  and  angle-of-sideslip  measurement  scale  factor 
Parameter  vector 
(PHI)  Roll  angle  (°) 
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Subscripts  (continued) 

a,  (5  Static  derivatives  with  respect  to  indicated  quantity 

<^a,  Control  derivatives  with  respect  to  indicated  quantity 

reference  Reference  aircraft  eg  condition:  longitudinal  -  45%c, 

vertical  -  177.5  in  above  datum  waterline,  lateral  -  on  the 
plane  of  symmetry 

UNITS 

It  is  a  common  convention  for  aircraft  operators  to  express  aircraft  weight,  altitude,  rate  of  ascent 
or  descent,  airspeed,  and  range  in  non  SI  units.  These  quantities  are  expressed  in  non  SI  units  in 
this  report  in  order  to  conform  with  this  convention  and  to  provide  uniformity  with  the  source 
data  on  which  the  results  are  based. 


Lift 


Positive  forces,  moments,  flow  angles  and  control  deflections  are  shown. 
Note:  Positive  control  deflections  give  negative  moments. 


vii 


Figure  1:  Force  and  moment  convention  (stability  axes) 
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1  INTRODUCTION 


The  General  Dynamics  F-lllC  aircraft  operated  by  the  Royal  Australian  Air  Force  (RAAF)  has 
a  different  geometry  from  the  versions  operated  by  the  United  States  Air  Force  (USAF).  The 
aerod5mamic  data  supplied  to  Australia  with  the  F-lllC  were  adapted  from  data  developed  for 
the  F-lllA  aircraft  which  has  the  same  fuselage,  and  for  the  F-lllB  which  has  the  same  wing 
[1],  [2].  To  validate  this  aerodynamic  data ,  a  flight  test  program  was  conducted  at  the  RAAF's 
Aircraft  Research  and  Development  Unit  (ARDU)  to  record  the  flight  behaviour  of  an  F-lllC 
subjected  to  manoeuvres  designed  to  optimise  the  identifiability  of  aerodynamic  characteristics. 
Analysis  of  the  data  has  been  carried  out  at  the  Air  Operations  Division  (AOD)  to  determine  the 
aerodynamic  characteristics.  Data  from  the  program  has  been  used  to  validate  a  comprehensive 
flight  dynamic  model  of  the  F-lllC  which  has  been  developed  at  AOD.  Data  from  the  model  is 
currently  being  used  in  support  of  the  F-lllC  Simulator  and  Avionics  Update  programs. 

With  the  introduction  of  the  variable  wing  sweep  F-lllC  into  service  in  1973,  the  RAAF  was 
faced  with  providing  technical  support  for  a  range  of  complex  systems.  In  particular,  the  air¬ 
craft  has  a  complex  adaptive  flight  control  system  [3]  which  is  required  to  provide  satisfactory 
flying  qualities  across  the  aircraft's  large  operational  flight  envelope.  At  the  time  of  acquisition 
it  was  not  possible  to  acquire  from  the  manufacturers  a  flight  dynamic  model  of  the  aircraft  suit¬ 
able  for  flight  dynamic  analysis.  To  provide  these  capabilities  a  model  was  developed  at  AOD 
from  the  manufacturer's  design  data  reports  and  from  wind  tunnel  tests  conducted  in  the  AOD 
Transonic  Wind  Tunnel.  The  model  allows  the  investigation  of  the  effects  of  pilot  commands, 
external  stores,  control  system  modifications  and  atmospheric  turbulence  on  the  aircraft's  rigid 
body  response  [4]. 

The  installation  of  flight  test  instrumentation  into  an  F-lllC  aircraft  at  ARDU  for  stores  clear¬ 
ance  testing  provided  an  opportunity  to  obtain  aerod5mamic  stability  and  control  data.  A  test 
program  was  developed  between  AOD  and  ARDU  and  was  approved  in  1983.  New  methods  of 
determimng  aerodynamic  data  from  flight  test  measurements  were  investigated  by  AOD  and 
sophisticated  software  for  this  purpose  was  acquired  from  the  National  Aeronautics  and  Space 
Admimstration  (NASA)  Dryden  Flight  Research  Center.  Additional  analysis  software  was  also 
developed  at  AOD.  Flight  testing  of  the  aircraft  was  carried  out  at  ARDU  in  1987  and  the  anal¬ 
ysis  of  the  data  was  conducted  at  AOD. 

This  report  is  one  of  a  series  of  F-lllC  aerodynamic  flight  data  analysis  reports,  and  presents  the 
longitudinal  and  lateral  stability  and  control  aerodynamic  derivatives  for  take-off  and  landing 
configurations.  Section  2  presents  details  of  the  flight  test  program.  Sections  3  and  4  describe 
the  flight  data  processing  and  data  analysis  procedures  and  section  5  discusses  the  results  of  the 
analysis.  Detailed  data  reduction  and  analysis  procedures  are  documented  in  [5]. 


2  FLIGHT  TEST  PROGRAM 


Testing  to  determine  the  stability  and  control  aerodynamic  derivatives  was  originally  plarmed 
under  ARDU  Test  Schedule  1667.  The  task  was  later  redefined  and  conducted  under  Test  Sched¬ 
ule  1691  to  take  into  account  variations  in  the  instrumentation  installation  compared  with  the 
installation  originally  planned.  Because  of  a  delay  in  the  manufacture  of  the  Nose  Boom  Trans- 
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Figure  2:  Instrumented  F-lllC  aircraft  (A8-132)  operated  by  ARDU  for  flight  trial 


ducing  Unit  (NBTU),  the  flight  testing  was  carried  out  in  three  phases.  Initial  flight  tests,  re¬ 
ferred  to  as  Phase  1,  were  conducted  during  February  and  March  1987  without  the  NBTU  fitted. 
Air  data  and  flow  direction  information  was  recorded  from  the  aircraft  Central  Air  Data  System 
(CADS).  The  NBTU  was  fitted  to  the  aircraft  in  September  1987.  Phase  2  consisted  of  a  number 
of  instrumentation  shakedown  flights.  Phase  3  testing  covered  additional  data  points  together 
with  repeated  test  points  from  Phase  1  and  these  were  carried  out  in  September  and  October 
1987. 

The  measurement  of  pressure  error  corrections  for  the  NBTU  was  carried  out  by  ARDU  as  part 
of  the  Test  Schedule  1691.  Details  of  these  measurements  together  with  other  details  of  the  test 
program  covering  the  aircraft,  test  equipment  and  the  scope  of  the  tests  are  presented  in  [6].  A 
summary  of  this  information  is  included  in  this  section. 


2.1  Test  Aircraft  and  Instrumentation 

The  F-lllC  test  aircraft  A8-132,  illustrated  in  figure  2,  was  extensively  modified  rmder  ARDU 
Test  Schedule  1650  with  the  addition  of  flight  test  quality  instrumentation  and  data  recording 
equipment.  This  equipment,  known  as  the  Airborne  Flight-Test  Recording  and  Analysis  Sys¬ 
tem  (AFTRAS),  provided  on-board  digital  magnetic-tape  recording  and  telemetry  information 
for  real-time  flight  test  monitoring.  The  instrumentation  was  developed  for  use  in  store-carriage 
and  release  tests  and  for  the  flight  dynamic  measurements.  Special  equipment  was  developed 
for  the  recording,  and  manual  adjustment  of  the  pitch  and  roll  adaptive  gain  settings  of  the  air¬ 
craft's  adaptive  control  system.  Due  to  a  lack  of  time,  instrumentation  required  to  monitor  other 
parts  of  the  adaptive  control  system  was  not  installed.  Similarly,  instrumentation  for  monitor¬ 
ing  engine  parameters  required  for  detailed  performance  measurements  was  not  fitted.  This 
instrumentation  was  not  required  for  the  estimation  of  the  aerodynamic  derivatives,  but  was 
required  for  the  related  investigation  of  the  adaptive  control  system  behaviour  and  aircraft  per¬ 
formance  characteristics.  As  a  result,  these  studies  were  not  completed. 

The  NBTU  used  for  Phase  3  of  the  trials  was  constructed  by  the  then  Advanced  Engineering  Lab¬ 
oratory  (AEL)  at  Salisbury,  and  was  designed  to  provide  pitot  pressure,  angle-of-attack,  angle- 
of-sideslip  and  linear  accelerations  parallel  and  normal  to  the  local  airflow  direction.  The  unit 
was  modelled  on  the  CONRAC®  Nose  Boom  Instrumentation  Unit  Model  25126F,  developed 
by  the  USAF  for  flight  dynamics  and  performance  measurements.  The  NBTU  is  shown  in  fig¬ 
ure  3  and  a  detailed  description  of  the  assembly  is  given  in  [7]. 
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Figure  3:  Nose  Boom  Transducing  Unit  (from  Reference  [7]) 


A  list  of  the  instrumentation  channels  (measurands),  including  their  range  and  resolution  [6], 
which  were  used  for  the  flight  dynamic  analysis  is  given  in  table  1.  Each  channel  was  recorded 
at  a  rate  of  60  samples  per  second. 


2.2  Weight  and  Centre  of  Gravity 

The  test  aircraft  was  weighed  prior  to  the  beginning  of  the  tests  and  had  a  zero  fuel  weight  of 
50  697  Ibf  with  the  eg  located  at  39.9%  mean  aerodynamic  chord  for  a  16°  wing  sweep  angle 
[6].  During  the  test,  the  aircraft  loading  and  eg  position  varied  in  accordance  with  the  aircraft 
auto-fuel  schedule.  The  forward  and  aft  fuel  contents  gauges  were  calibrated  prior  to  the  test 
program  and  were  manually  recorded  at  each  test  point.  During  the  test  manoeuvres,  the  air¬ 
craft  weight  was  typically  in  the  region  of  70  000  Ibf. 


2.3  Aircraft  Test  Configurations 

In  this  report,  results  are  presented  for  the  aircraft  in  take-off  and  landing  configurations  for 
nominal  16°  and  26°  wing  sweep  angles.  The  take-off  and  landing  configurations  consisted  of: 

•  landing  gear  down, 

•  wing-root  glove  rotated  (trailing  edge  up). 
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Quantity 

Symbol 

Units 

Range 

Resolution 

ACCEL  LONG  CG* 

g 

±5 

±0.05 

ACCEL  LAT  CG 

ay 

g 

±5 

±0.05 

ACCEL  VERT  CG*^ 

0-z 

g 

±10 

±0.05 

ROLL  RATE 

P 

% 

±300 

±3 

PITCH  RATE 

q 

% 

±100 

±1 

YAW  RATE 

r 

7s 

±50 

±0.5 

ROLL  ACC  CG 

P 

rad/s^ 

±10 

±0.05 

PITCH  ACC  CG 

4 

rad/s^ 

±5 

±0.05 

YAW  ACC  CG 

r 

rad/s^ 

±5 

±0.05 

ROLL  ANGLE 

O 

±180 

±0.5 

PITCH  ANGLE 

9 

0 

±180 

±0.5 

ANGLE-OF-ATTACK* 

a 

o 

-3  ^  25 

±0.5 

ANGLE-OF-SIDESLIP* 

P 

0 

±24 

±0.5 

VELOCITY 

V 

kn 

0  ^  900 

±10 

MACH  Nod 

M 

— 

0.3  1.3 

±0.001 

altitude! 

H 

ft 

-500  -4  5  5000 

±1.5 

WING  SWEEP 

A 

0 

16  72.5 

±0.05 

STABILATOR  (right) 

^StR 

o 

-30  ^  15 

±0.1 

STABILATOR  (left) 

o 

-30  ->  15 

±0.1 

RUDDER 

Sr 

o 

±30 

±0.1 

SPOILER  (right) 

o 

0  ^  45 

±0.1 

SPOILER  (left) 

4spr 

0 

0  45 

±0.1 

STICK  POS  (longitudinal) 

in 

-4.4  ^  3.6 

±0.05 

STICK  POS  (lateral) 

in 

±5 

±0.05 

RUDDER  PED  POS 

in 

±3 

±0.03 

*  NBTU  measurements  were  also  available  in  Phase  3 
t  fine  and  coarse  readings  available 

I  crew  module  normal  accelerometer  used  in  Phase  3  for  the 
later  part  of  flight  1  and  all  of  flights  2  and  3  because  of  a 
fault  with  the  eg  accelerometer  signal 

Table  1:  Instrumentation  channels  used  for  flight  dynamic  analysis 
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•  leading-edge  slats  extended  and 

•  flaps  deflected  (0°,  15°,  25°  and  35°). 

The  configurations  are  summarised  in  table  2. 


2.4  Flight  Control  System  Status 

The  tests  were  conducted  with  the  flight  control  system  in  its  normal  mode  and  with  the  system 
gains  determined  by  the  adaptive  mode  gain  changer.  However,  when  some  manoeuvres  were 
performed,  motion  due  to  the  adaptive  mode  was  superimposed  onto  the  natural  motion  of  the 
aircraft.  In  these  cases,  the  contribution  from  the  adaptive  mode  was  decreased  by  the  pilot  by 
pumping  the  stick  to  reduce  the  gain. 


2.5  Test  Points 

The  matrix  of  take-off  and  landing  configuration  test  points  covered  in  Phase  1  and  Phase  3  is 
given  in  table  2.  It  should  be  noted  that  the  take-off  and  landing  configuration  tests  were  per¬ 
formed  at  a  nominal  altitude  of  1000  ft  and  were  thus  out  of  ground  effect. 

In  Phase  3,  a  number  of  Phase  1  test  point  manoeuvres  were  repeated  to  compare  the  accuracy  of 
CADS  measurements  relative  to  the  NBTU  system  measurements  of  angle-of-attack  and  angle- 
of-sideslip. 


2.6  Test  Manoeuvres 

The  following  manoeuvres  were  performed  at  each  test  point: 

1.  Accurate  trim. 

2.  Pitch  input  (from  1  g  to  approximately  2  g  and  return  to  1  g). 

3.  Trim. 

4.  Pitch  input  (from  1  g  to  approximately  0  g  and  return  to  1  g). 

5.  Accurate  trim. 

6.  Rudder  step  input  to  left  followed  by  aileron  doublet  to  achieve  =f30°  bank  angle.  Rud¬ 
der  and  aileron  released  together. 

7.  Trim. 

8.  Manoeuvre  6  repeated,  but  with  rudder  input  to  right  and  opposite  roll  applied. 

9.  Trim. 
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Manoeuvre 

Phase 

Speed 

Sweep 

Manoeuvre  Name 

(KCAS) 

(°) 

Longitudinal  Lateral 

Landing  Configuration 


15°  flap 

1 

177 

16 

plf6el0 

plf6el4 

188 

23 

plf6e26 

plf6e30 

236 

16 

plf6e01 

plf6e05 

238 

23 

plf6el8 

plf6e22 

34°  flap 

1 

147 

16 

plf6e50 

plf6e54 

150 

20 

plf6e58 

plf6e62 

3 

200 

16 

p3fle72 

p3fle74 

160 

26 

p3fle78 

p3fle80 

Slats,  no  flap 

1 

198 

16 

plf6e42 

plf6e46 

215 

25 

plf6e34 

plf6e38 

240 

16 

plf7e77 

plf7e81 

240 

26 

plf7e69 

plf7e73 

Take-off  Configuration 


25°  flap 

1 

158 

16 

plf7e85 

plf7e89 

3 

165 

16 

p3fle84 

p3fle86 

Table  2:  Take-off  and  landing  test  point  matrix  for  Phases  1  and  3,  altitude  =  1000  ft 


Each  manoeuvre  is  referred  to  by  a  simple  naming  convention  as  shown  in  table  2.  Each  ma¬ 
noeuvre  has  a  name  of  the  form  pifjekk[x],  where  i  indicates  the  Phase  of  the  tests,  j  the  flight 
number  and  kk  the  event  number  as  recorded  on  the  flight  tape.  An  optional  x  indicates  an 
'a',  'b'  or  'c'  manoeuvre.  The  first  longitudinal  and  lateral  manoeuvre  is  referred  to  as  the  'a' 
manoeuvre  and  the  second  the  ‘h'  manoeuvre.  A  'c'  manoeuvre  is  a  combined  'a'  and  'b  ma¬ 
noeuvre. 

One  noticeable  difference  between  the  behaviour  of  the  aircraft  in  the  clean  configuration  and 
the  take-off  and  landing  configuration  was  the  average  corrected  sideslip  excursion.  In  the  clean 
configuration  this  was  approximately  ±6°,  while  in  the  take-off  and  landing  configuration  it  in¬ 
creased  to  approximately  ±13°.  Typical  longitudinal  and  lateral  aircraft  inputs  and  responses 
are  shown  in  figures  4, 5, 6  and  7.  The  effect  of  the  flight  control  system  is  pronounced,  as  illus¬ 
trated  in  figures  4  and  6,  where  comparisons  of  stick  and  pedal  position  are  made  with  their 
respective  control  surface  deflections. 

3  FLIGHT  DATA  PROCESSING 
3.1  AFTRAS  Data  Extraction 

The  AFTRAS  flight  data  system  provided  procedures  for  accessing  selected  measurands,  for  ap¬ 
plying  calibrations  to  give  engineering  units,  and  for  formatting  the  data  for  subsequent  anal¬ 
ysis.  Calibrations  were  stored  as  polynomial  coefficients  and  included  the  date  of  calibrations 
and  were  appended  to  the  data  files  obtained  from  each  flight.  Details  of  these  procedures  are 
given  in  [8].  The  procedures  were  carried  out  using  program  EXTRACT  which  was  run  on  DEC 
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Figure  4:  Longitudinal  control  input  and  aircraft  response  for  case  plf6e01a 
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Figure  5:  Longitudinal  aircraft  response  for  case  plf6e01a 


DSTO-TR-0321 


Wing  sweep  =  26  deg.  Altitude  =  1300  ft  223  KCAS  0  deg.  flap 


Fri  Sep  23  16:49:16  1994 


-10.00  -i  . . . 

-20,00  -j  .  ; .  i . 1 

'  '  I  , 

0.00  5.00  10.00  15  00 

TIME.1(00:06:55.0092)  (s) 


30.00  • 
20.00  ■ 
10.00  - 
0.00- 
-10.00  - 
-20.00  - 
-30.00  - 

0.00 


. / 

/ 

\^' . 

_ _ _ 

. / 

. \.. 

v/ 

i 

■(|)(o) 


T1ME.1(00:06:55.0092) 


(S) 


1  ■  -/ . 

. 

j 

^  :  1  / 

1/ 

. 

. 

'll 

0.00  5.00  10.00  15.00 

TIME.1(00:06:S5,0092)  TIME  (s) 


6r(“) 

RPED*10  (in) 


LATSTK-IO(in) 


Figure  6:  Lateral  control  input  and  aircraft  response  for  case  plf6e38a 
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VAX®  computers  at  ARDU  or  the  Aeronautical  and  Maritime  Research  Laboratory  (AMRL). 

3.2  Data  Corrections 

The  techniques  for  determining  aerodynamic  parameters  from  the  flight  time  histories  are  de¬ 
scribed  in  section  4.  The  data  for  input  to  the  parameter  estimation  software  needs  additional 
processing  for  conversion  to  the  required  format  and  to  apply  measurement  error  corrections. 
This  processing  stage  was  carried  out  at  AOD  with  the  program  PDF  on  an  IBM  UNIX  ®  work¬ 
station  and  is  detailed  in  [5].  The  operations  performed  were  : 

1.  to  correct  the  static  pressure  measurement  for  sensor  position  error  [6],  [5], 

2.  to  adjust  the  airspeed  and  Mach  number  for  scale  altitude  [5], 

3.  to  apply  time  shifts  to  the  time  history  records  to  allow  for  instrument,  signal  condi¬ 
tioning  and  recording  lags  (see  section  3.3), 

4.  to  adjust  the  angle-of-attack  and  angle-of-sideslip  measurements  for  sensor  location 
(see  section  3.4), 

5.  to  calculate  elevator  and  aileron  deflections  from  the  measured  stabilator  deflectionis, 
and 

6.  to  calculate  weight,  eg  and  inertia  information  from  the  fuel-used  data. 

3.3  Instrumentation  Time  Lags 

A  procedure  for  identifying  the  relative  time  lags  between  instrumentation  channels  was  de¬ 
veloped  at  AOD  and  is  documented  in  [9].  This  procedure  uses  the  maximum  likelihood  tech¬ 
nique  and  was  applied  to  a  number  of  selected  time  histories  to  determine  the  lag  parameters. 
The  resulting  parameters  are  given  in  table  3  where  a  positive  integer  indicates  n  time  samples 
(equivalent  to  a  time  interval  of  n/60  seconds)  lagged  with  respect  to  the  control  deflection  sig¬ 
nals.  The  table  shows  that  the  NBTU  signal  for  angle-of-attack  in  fact  leads  the  control  deflection 
signals  by  two  sample  intervals.  The  time  lag  differences,  for  olcads  and  6,  between  the  Phase  1 
and  Phase  3  flights  are  most  likely  due  to  changes  to  the  instrumentation  during  the  installation 
of  the  NBTU. 

3.4  Calibration  of  a.  and  /3  Flow  Vanes 

It  is  necessary  to  determine  the  position  errors  and  calibration  constants  for  the  angle-of-attack 
and  angle-of-sideslip  transducers.  In  particular  the  aircraft  CADS  transducers  were  mounted 
close  to  the  forward  fuselage  where  local  flow  angles  can  differ  substantially  from  free-stream 
values.  A  method  for  determining  these  calibration  constants  is  described  in  [10].  This  method 
uses  combined  parameter  and  state  estimation  techniques.  The  estimation  software  which  is  de¬ 
scribed  in  [11]  and  [12]  was  developed  for  AOD  imder  a  research  agreement  with  the  University 
of  Newcastle. 
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Channel 

Phase  1 

Phase  3 

0 

0 

0 

0 

Sr 

0 

0 

SspR 

0 

0 

Sspi, 

0 

0 

OtCADS 

2 

-2 

OiNBTU 

na 

-2 

P 

0 

0 

9 

0 

-2 

0 

0 

V 

2 

2 

Q. 

3 

3 

r 

2 

2 

q 

na 

3 

an 

4 

4 

nu 

0 

^  ^  na  -  not  available 
Note  ,  , 

nu  -  not  used 


Table  3:  Instrumentation  lags 


Application  of  the  method  to  the  measured  time  histories  gave  estimates  for  the  calibration  con¬ 
stants  for  the  CADS  and  NBTU  angle-of-attack  and  sideslip  measurement  systems.  The  calcula¬ 
tions  of  the  angle-of-attack  and  sideslip  vane  calibration  constants  are  described  in  appendices 
A  and  B  respectively. 

For  the  Phase  1  angle-of-attack  sensor  it  was  found  that  the  calibration  constant  was  invariant 
with  Mach  number.  A  constant  value  of  1.04  for  Phase  1  was  used.  During  Phase  3  the  limited 
number  of  test  cases  showed  that  the  calibration  constant  was  not  invariant  with  Mach  number 
and  each  individual  value  was  used  in  the  estimation  procedure. 

Similar  behaviour  was  observed  for  the  angle-of-sideslip  sensor,  with  a  constant  value  of  1.55 
being  used  for  Phase  1  and  the  individual  values  for  each  of  the  three  Phase  3  flights.  The  NBTU 
calibration  constants  are,  as  expected,  smaller  than  the  CADS  values  because  the  NBTU  is  sens¬ 
ing  the  sideslip  flow  angles  in  a  relatively  undisturbed  flow,  compared  to  the  CADS  sensor  on 
the  imderside  of  the  fuselage. 

The  angle-of-attack  and  sideslip  signals  must  be  corrected  also  because  the  vanes  are  not  located 
at  the  aircraft  eg  and,  therefore,  aircraft  rotations  are  incorrectly  sensed  by  the  vanes  as  angle-of- 
attack  or  sideslip  changes.  The  corrections  are  given  in  appendices  A  and  B.  It  should  be  noted 
that  the  angle-of-attack  and  sideslip  plots  in  figures  4,  6,  8  and  13  are  for  the  raw  measured 
signal,  not  the  corrected  signal. 
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4  DATA  ANALYSIS  PROCEDURES 

4.1  Aerodynamic  Parameter  Estimation 


Several  techniques  have  been  developed  in  the  last  two  decades  for  the  estimation  of  aerody¬ 
namic  derivatives  from  flight  test  measurements.  The  aerod3mamic  and  control  derivatives  of 
the  F-lllC  aircraft  were  estimated  using  the  output-error  parameter  estimation  technique  de¬ 
veloped  at  NASA  Dryden  Flight  Research  Facility  [13].  This  technique  involves  the  calculation 
and  minimisation  of  a  cost  function.  The  cost  function  is  the  weighted  sum  of  the  square  of 
the  difference  between  the  computed  model  response  and  the  flight  measured  response.  The 
technique  and  associated  computer  program,  pEst,  are  described  fully  in  [14]  and  briefly  in  the 
following  sections. 

It  should  be  noted  that  in  [15],  [16]  and  [17],  a  maximum  likelihood  implementation  of  the 
output-error  technique  was  used.  In  this  implementation,  the  response  weighting  matrix  is  ad¬ 
justed  to  provide  an  optimally  minimal  cost  function.  The  associated  software,  MMLE3  [18]  has 
been  replaced  by  pEst  which  has  fixed  weightings  on  the  response  variables.  Additionally,  pEst 
provided  a  far  more  user-friendly  interface,  allowing  on-line  interaction  during  the  analysis  as 
is  shown  in  appendices  C  and  D  where  example  pEst  runs  are  given  for  both  a  longitudinal  and 
a  lateral  test  case. 


4.2  Equations  of  Motion 

Within  the  range  of  flight  conditions  tested  in  the  validation  of  the  F-lllC  flight  dynamic  model, 
it  has  been  assumed  that  the  aircraft  motion  can  be  adequately  represented  by  separate  classi¬ 
cal  flight  dynamic  models  for  longitudinal  and  lateral  motion.  The  manoeuvres  have  been  de¬ 
signed  so  that  the  longitudinal  and  lateral  motion  of  the  aircraft  may  be  analysed  separately 
to  give  a  set  of  stability  and  control  derivatives  describing  the  longitudinal  and  lateral  aerody¬ 
namic  forces  and  moments. 

The  flight  dynamic  model  used  for  this  purpose  was  the  standard  model  described  in  [14]  which 
assumes  small  disturbance  motion  and  linear  aerodynamic  characteristics.  It  incorporates  the 
non-linear  equations  of  motion  presented  in  appendices  E  and  F. 


4.3  A  Priori  Data 

A  priori  information  can  be  used  to  assist  the  identification  process.  Previous  estimates  of  aero¬ 
dynamic  derivatives  are  often  available  from  sources  such  as  wind  tuimel  data  reports.  Param¬ 
eter  estimation  results  from  previous  similar  cases  may  also  be  used  to  speed  up  convergence. 
In  cases  where  a  derivative  makes  only  a  small  contribution  to  the  motion,  parameters  may  be 
constrained  to  their  a  priori  values.  The  a  priori  values  were  obtained  from  the  wind  tunnel 
values  used  in  the  imvalidated  flight  dynamic  model  [4]  discussed  in  section  1. 
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4.4  Cramer-Rao  Bounds 

For  the  unconstrained  parameters  the  estimation  procedure  calculates  a  measure  of  the  estima¬ 
tion  accuracy  known  as  the  Cramer-Rao  bound.  The  interpretation  of  this  quantity  is  given  in 
[18].  The  results  are  plotted  for  each  aerodynamic  derivative  and  can  be  used  in  conjunction 
with  the  observed  repeatability  to  indicate  the  estimation  accuracy.  To  account  for  the  fact  that 
the  measurement  signal  noise  is  bandwidth  limited  compared  to  the  theoretical  assumption  of 
white  noise,  i.e.  infinite  bandwidth,  the  Cramer-Rao  bound  is  factored  by  a  multiple  of  10  in 
accordance  with  the  procedures  described  in  [19]. 


4.5  Cost  Function 

The  cost  function  is  an  explicit  function  of  the  computed  response,  and  thus  an  implicit  function 
of  the  vector  of  unknown  parameters.  The  cost  function  used  in  the  program  is  : 


1  rit 

i=:l 

where  z  is  the  response  computed  by  integrating  the  equations  of  motion,  z  is  the  flight  mea¬ 
sured  response,  ^  is  the  parameter  vector,  ut  and  are  the  number  of  time  history  and  response 
variables  respectively.  FF  is  a  response  weighting  matrix  which  is  defined  by  the  user  and  gives 
an  indication  of  confidence  in  each  measurement. 

The  matches  between  the  computed  responses  and  the  flight  measured  responses  in  the  take-off 
and  landing  configuration,  obtained  with  the  default  weightings,  and  those  used  for  the  clean 
aircraft  configuration  results,  [15],  [16],  [17],  [20],  [21],  [22],  were  poor.  The  aerodynamic 
derivatives  estimated  also  compared  poorly  with  the  expected  results.  Some  of  the  damping 
derivatives  had  the  opposite  sign  than  what  was  expected,  incorrectly  indicating  that  the  air¬ 
craft  was  unstable  in  certain  modes.  Other  derivatives  possessed  unreasonably  high  magni¬ 
tudes.  The  response  weightings  are  a  function  of  the  range  of  the  response  measured.  Upon 
investigation  it  was  found  that  the  take-off  and  landing  manoeuvres  showed  an  increase  in  the 
range  of  some  of  the  measurements  compared  with  the  clean  aircraft  case,  e.g.  angle-of-sideshp. 

A  scheme  for  establishing  a  new  set  of  response  weightings  was  therefore  devised.  The  scheme 
made  use  of  the  full  range  of  the  measured  response  and  relative  confidence  in  each  measure¬ 
ment  response.  The  half  range  of  each  measured  response  was  inverted  and  multiplied  by  100 
to  give  the  initial  matrix  of  weights.  The  factor  of  100  was  used  solely  to  obtain  a  matrix  of  whole 
numbers.  Then  a  confidence  factor,  ranging  from  1.0  to  2.0,  was  assigned  to  each  measured  re¬ 
sponse.  For  example  a  value  of  2.0  indicated  high  relative  confidence  in  the  measured  response. 
The  initial  response  weighting  matrix  was  then  multiplied  by  the  matrix  of  confidence  factors 
to  obtain  the  final  weighting  matrix.  Tables  4  and  5  summarise  the  procedure  for  establishing 
the  longitudinal  and  lateral  weighting  matrices  respectively.  It  should  be  noted  that  the  lateral 
response  weighting  matrix  required  further  refinement  after  preliminary  pEst  runs  indicated 
less  than  satisfactory  matches  between  the  measured  and  calculated  responses.  The  sensitivity 
of  the  resultant  matches  to  weighting  matrix  variations  is  indicative  of  the  tow  confidence  in 
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Response 

Range 

Inverse'^100 

Confidence 

Weight 

a 

Oto  11 

18.20 

1.5 

27 

q 

±15 

6.67 

2.0 

13 

e 

Oto  15 

13.00 

1.0 

13 

On 

0  to  2 

100.00 

2.0 

200 

Table  4:  Longitudinal  response  weighting  matrix. 


Response 

Range 

Inverse*100 

Confidence 

Weight 

Adjusted 

Confidence 

Adjusted 

Weight 

P 

±10 

10.0 

1.5 

15.0 

0.75 

7.5 

P 

±40 

2.5 

2.0 

5.0 

1.76 

4.4 

r 

±8 

12.5 

2.0 

25.0 

3.20 

40 

</> 

±40 

2.5 

1.0 

2.5 

1.60 

4 

Oy 

±0.3 

333.0 

2.0 

666.0 

4.00 

1333 

Table  5:  Lateral  response  weighting  matrix. 


some  of  the  estimated  parameters  primarily  because  of  the  flight  control  system  damping  out 
the  natural  response  of  the  aircraft.  It  should  be  noted  that  the  estimated  derivatives  themselves 
are  insensitive  to  changes  in  the  overall  weightings.  Experience  at  NASA  Dryden  has  indicated 
that  doubling  of  the  weighting  should  result  in  a  5%  or  less  change  in  parameters.  The  purpose 
of  calculating  a  new  set  of  weightings  was  to  obtain  a  better  match  between  the  calculated  air¬ 
craft  responses  and  the  flight  measured  responses.  The  adjusted  weighting  matrix  is  shown  in 
table  5. 


4.6  Reference  eg  Correction 


Longitudinal,  and  to  a  lesser  extent,  lateral  derivatives,  vary  with  eg  position.  The  eg  position 
varies  with  fuel  usage  and  the  flight  test  results  were  calculated  at  the  actual  flight  test  eg  posi¬ 
tion.  The  wind  tunnel  results  have  been  plotted  with  respect  to  the  flight  test  eg  as  well  as  the 
reference  eg  position.  The  reference  eg  was  at  45%  mean  aerodynamic  chord,  177.5  in  above  the 
waterline  and  at  the  zero  reference  point  laterally. 

A  full  summary  of  the  eg  corrections  is  given  in  [23]. 
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5  RESULTS  AND  DISCUSSION 
5.1  Longitudinal  Derivatives 

For  each  longitudinal  manoeuvre  the  output-error  parameter  estimation  program,  pEst,  was 
used  with  the  elevator  deflection  as  input,  and  with  the  following  measurands  as  responses  : 

•  angle-of-attack,  a, 

•  pitch  rate,  q, 

•  pitch  attitude,  6,  and 

•  normal  acceleration,  a^. 

In  the  clean  configuration  analysis  [15],  [16],  [17],  the  pitch  attitude  was  recorded  but  was  given 
zero  weighting  in  the  identification  process.  For  the  take-off  and  landing  configurations,  the 
pitch  rate  data  alone  did  not  provide  adequate  information  for  the  identification  process.  Hence 
the  pitch  attitude  was  also  included  as  a  response.  In  general,  for  all  derivative  identification  it 
was  found  that  the  best  results  were  obtained  when  the  aircraft  was  trimmed  accurately  before 
each  manoeuvre.  As  was  shown  in  [15],  [16]  and  [17],  on  some  occasions  trim  was  not  ob¬ 
tained  between  consecutive  longitudinal  manoeuvres,  so  for  these  cases  the  two  manoeuvres 
were  analysed  as  one.  This  resulted  in  an  overall  reduction  in  the  number  of  test  cases.  Allow¬ 
ing  the  aircraft  to  trim  after  each  manoeuvre  enabled  an  accurate  'end  of  manoeuvre'  time  to  be 
defined.  Analysing  the  data  beyond  this  point  resulted  in  a  reduction  in  the  consistency  of  the 
derivatives. 

An  example  of  the  time  histories  is  given  in  figure  8.  The  flight  measured  response  is  drawn 
as  a  solid  line.  The  computed  response  is  drawn  as  a  dotted  line  and  is  also  identified  by  the 
addition  of  -HAT  to  the  name.  The  time  history  plots  in  figure  8  and  the  Cramer-Rao  boimds 
were  used  to  assess  the  accuracy  of  the  identification  procedure. 

Initially,  the  responses  were  not  well  matched  for  a  number  of  cases.  The  matches  were  im¬ 
proved  by  adjusting  the  response  weighting  as  discussed  in  section  4.5.  The  normal  accelera¬ 
tion  weighting  was  increased  significantly  from  the  default  pEst  value,  improving  the  match, 
and  this  led  to  a  closer  match  between  the  flight  estimated  normal-force-coefficient  derivative 
values  and  the  wind  tunnel  results. 

The  pEst  program  calculates  the  coefficients  in  body  axes.  To  convert  to  stability-axis  coefficients 
the  following  relationship  is  used  : 


Cl  =  Cncosol  —  CAsina 

It  should  be  noted  that  the  thrust  was  not  measured  in  the  flight  test  program.  Thus  none  of 
the  axial  force  derivatives  have  been  estimated.  The  error  introduced  by  neglecting  the  axial 
force  derivatives  is  small  for  low  angles-of-attack.  See  [24]  for  a  more  detailed  description  of 
stability-axis  and  body-axis  systems. 

The  longitudinal  derivatives  identified  are  given  in  table  6  and  the  results  shown  in  figures  9, 
10,  11  and  12.  The  derivatives  estimated  from  the  flight  tests  are  for  the  actual  flight  eg.  To 
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Normal  Force 

Pitching  Moment 

Aerodynamic 

C'aTc 

Crr. 

Cr„ 

Control 

Cn.. 

^rns. 

Table  6:  Longitudinal  derivatives  estimated 


enable  valid  comparisons  the  wind  tunnel  data  have  been  adjusted  to  the  flight  eg  value.  The 
C'ATi/  and  Cm^  aerodynamic  derivatives  used  in  the  longitudinal  flight  dynamic  model  were 
not  identified  due  to  their  small  contribution  and/ or  the  difficulty  in  identification. 


5.1.1  Angle-of-Attack  Derivatives 

GenerallytheC'i\^„  derivatives  shown  in  figures  9,  10,  Hand  12  display  reasonable  agreement 
with  wind  tunnel  results  for  the  same  flight  conditions,  including  flap  and  eg  position.  How¬ 
ever,  the  value  estimated  for  the  26°  wing  sweep  case  with  34°  flap  angle,  shown  in  figure 
12,  appears  to  be  poor  and  the  large  Cramer-Rao  bound  reflects  this. 

The  pitch  stiffness  derivative,  Cma>  is  a  linear  fxmetion  of  eg  position.  Figures  11  and  12  show 
that  the  estimated  Cm^  values,  for  both  wing  sweep  angles,  are  significantly  lower  in  magnitude 
than  the  wind  tunnel  values  adjusted  to  the  flight  eg  and  that  the  pitch  stiffness  increases  with 
flap  deflection.  It  should  be  noted  that  for  the  16°  wing  sweep  case  the  wind  tunnel  values  of 
pitch  stiffness  at  the  reference  eg  position  are  unstable.  This  is  because  the  reference  eg  is  an 
unrealistic  value  for  a  16°  wing  sweep;  correcting  the  eg  for  actual  flight  conditions  results,  as 
expected,  in  a  stable  derivative. 


5.1.2  Pitch  Rate  Derivatives 

During  the  initial  estimation  of  the  longitudinal  derivatives,  showed  considerable  variabil¬ 

ity  when  identified.  The  Cat,  derivatives  estimated  for  the  clean  aircraft  configuration  [15],  [16], 
[17]  were  in  good  agreement  with  the  wind  tunnel  data  and  so  wind  tunnel  values  of  Cat,  for 
the  take-off  and  landing  configuration  were  used  as  a  priori  data  and  Cat,  was  fixed.  Fixing  Cat, 
did  not  affect  the  final  match  nor  did  it  cause  any  significant  changes  in  the  remaining  estimated 
derivatives. 

The  total  restoring  pitching  moment  following  a  disturbance  consists  of  contributions  due  to 
the  stabilator  commanded  by  the  automatic  control  system  as  well  as  contributions  due  to  pitch 
rate,  angle-of-attack  and  rate  of  change  of  angle-of-attack.  An  ideal  manoeuvre  to  identify  the 
pitching  moment  derivatives  would  be  a  sudden  stabilator  (pulse  type)  input  followed  by  zero 
input  to  allow  the  aircraft  to  damp  naturally.  Unfortunately  in  some  manoeuvres  the  stabilator 
inputs  were  varying  over  the  entire  manoeuvre,  due  either  to  the  pilot  adjusting  the  stabilator 
position,  or  the  aircraft  control  system  providing  artificial  pitch  damping.  Although  the  identi¬ 
fication  model  is  still  applicable  in  these  cases,  the  resulting  response  contains  little  information 
from  which  to  determine  the  pitch  damping  derivatives  accurately.  Testing  with  the  automatic 
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Figure  10:  Longitudinal  aerodynamic  derivatives  for  A  =  26' 
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pitch  damping  control  system  switched  off  would  not  be  safe  because  of  the  resulting  deteri¬ 
oration  in  the  aircraft  handling  behaviour.  As  a  result  of  the  difficulties  involved  in  obtaining 
an  accurate  estimate  of  the  Cm^  derivative,  the  final  estimates  for  both  16°  and  26°  wing  sweep 
exhibited  a  greater  degree  of  variability  in  their  identification  than  the  other  longitudinal  deriva¬ 
tives. 

Crm,  and  (7^,  are  difficult  to  identify  separately.  However  it  should  be  noted  that  some  success 
has  been  met  using  the  maximum  likelihood  technique  with  specially  designed  manoeuvres, 
such  as  two-point  hesitation  rolls  [25].  In  [15],  [16]  and  [17],  the  and  derivatives 
were  estimated  using  the  MMLE3  'hard  constraints'  facility,  i.e.  the  value  of  one  derivative  was 
constrained  to  be  a  multiple  of  the  other.  was  estimated  from  Cm,  using  this  technique 
and  the  results  matched  the  wind  tunnel  data  well.  In  pEst  however,  there  is  no  provision  for 
such  'hard  constraints',  and  so  for  the  take-off  and  landing  configuration  cases  it  was  decided 
to  accept  the  wind  tunnel  values  and  identify  only  Cm,-  Earlier  estimation  carried  out  on 
the  clean  aircraft  configuration  [15],  [16],  [17],  had  shown  sufficient  agreement  between  the 
estimated  Cm^,  values  and  the  wind  tunnel  data  to  support  this  decision. 


5.1.3  Control  Derivatives 

In  a  similar  fashion  to  Cm,  and  the  derivatives  C^^  and  C^g^  were  estimated  for  the  clean 
configuration  [15],  [16],  [17],  using  the  MMLE3 'hard  constraint' facility  due  to  the  difficulty  in 
identifying  them  separately.  Again,  since  there  is  no  provision  for  'hard  constraints'  in  pEst,  Cj^^ 
was  fixed  using  the  wind  tunnel  values  as  a  priori  data.  The  measurement  of  dynamic  deriva¬ 
tives  in  wind  tunnels,  for  example  C^,,  is  generally  difficult  in  comparison  to  static  derivatives, 
such  as  but  as  with  Cm^  there  had  been  sufficient  agreement  between  the  estimated  Cm, 
values  and  the  wind  tunnel  data  for  the  clean  configuration  [15],  [16],  [17]  to  justify  this  deci¬ 
sion. 

For  the  16°  wing  sweep  case,  the  flight  test  results  for  Cms^  and  Cmg^  compare  well  with  the 
wind  tunnel  data,  although  the  Cramer-Rao  bounds  indicate  a  higher  level  of  confidence  in  the 
estimates  than  the  Cn{^  estimates.  At  26°  wing  sweep  the  flight  test  results  for  Cmg^  also 
compare  well  with  the  wind  tunnel  data,  differing  by  less  than  15%. 


5.2  Lateral  Derivatives 

The  output-error  parameter  estimation  program,  pEst,  was  applied  to  each  lateral  manoeuvre. 
For  wing  sweeps  less  than  47°  spoilers  are  automatically  deployed  for  roll  control  in  addition 
to  the  aileron  deflection  [4].  Hence  in  all  take-off  and  landing  configurations,  where  the  wing 
sweep  is  always  less  than  47°,  both  spoilers  and  ailerons  are  used  as  inputs.  The  following  mea- 
surands  were  used  as  responses  : 

•  angle-of-sideslip,  (3, 

•  roll  rate,  p, 

•  yaw  rate,  r. 
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Side  Force 

Rolling  Moment 

Yawing  Moment 

Aerodynamic 

Cvp 

Ch 

Cn, 

Cl. 

Cn, 

Control 

Cls. 

Cns^ 

Cls., 

Cns,f, 

Cls. 

Table  7:  Lateral  derivatives  estimated 

•  lateral  acceleration,  ay,  and 

•  roll  attitude,  4>. 

The  roll  attitude  {(p)  was  recorded  but  was  initially  given  zero  weighting  in  the  identification 
process  for  the  clean  aircraft  configuration  [20],  [21],  [22].  However,  acceptable  time  history 
matches  were  difficult  to  obtain  in  the  take-off  and  landing  configuration  without  cp  matching. 
The  match  improved  considerably  when  (p  was  included  as  an  output  match. 

An  example  of  the  time  histories  for  these  output  variables  for  a  given  rudder  and  aileron  input 
is  shown  in  figure  13.  The  time  history  plots  shown  in  figure  13  and  the  Cramer-Rao  bounds 
were  used  to  assess  the  accuracy  of  the  identification  procedure. 

The  lateral  derivatives  identified  are  given  in  table  7  and  the  results  shown  in  figures  14,  15, 
16  and  17. 


5.2.1  Angle-of-Sideslip  Derivatives 

The  static  lateral  aerod)mamic  derivatives  are  the  side  force,  rolling  and  yawing  moments  due 
to  sideslip,  and  Cn^  respectively.  For  the  side  force  due  to  sideslip  derivative  Cy^,  the 

flight  test  results  show  little  scatter  and  agree  well  with  the  wind  tunnel  results. 

The  derivative  Ci^,  known  as  the  dihedral  effect,  matches  the  wind  tunnel  data  well  for  16°  wing 
sweep.  At  26°  wing  sweep  the  magnitudes  of  the  flight  data  are  consistently  less  than  that  of 
the  wind  tunnel  values  by  approximately  30%.  The  Ci^  estimates  displayed  similar  behaviour 
in  [20],  [21]  and  [22]  for  the  clean  aircraft  configuration. 

The  weathercock  stability  derivative,  C„^,  results  shows  little  scatter  at  most  flight  conditions. 
For  16°  wing  sweep  the  flight  test  data  tends  to  be  about  40%  of  the  wind  tunnel  value.  At  26° 
wing  sweep  the  flight  data  tends  to  be  about  50%  of  the  wind  tunnel  value.  This  is  consistent 
with  the  results  given  in  [20],  [21]  and  [22]  where  the  wind  tunnel  results  regularly  overesti¬ 
mated  values  of  Cn^- 


5.2.2  Roll  and  Yaw  Rate  Derivatives 

The  aerod5mamic  derivatives  with  the  largest  identification  variability  are  the  side  force  due  to 
roll  and  yaw  rate  derivatives,  Cy^,  and  Cy;,  the  rolling  moment  due  to  yaw  rate  derivative  Ci,., 
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and  the  yawing  moment  due  to  roll  rate  derivative  .  In  general  for  conventional  aircraft  such 
as  the  F-lllC,  and  Cy^  are  not  physically  significant,  whereas  the  cross  derivatives  Ci^  and 
Cnp  may  be  significant  but  the  responses  from  the  lateral  manoeuvres  contain  little  information 
to  determine  these  derivatives.  Therefore  they  were  not  estimated  but  were  fixed  at  the  a  priori 
values. 

Although  the  match  obtained  using  the  a  priori  data  was  quite  acceptable,  some  of  the  identified 
parameters  were  obviously  incorrect.  For  example,  the  estimated  value  of  the  non-dimensional 
yawing  moment  due  to  yaw  rate  derivative  Cn^,  which  is  a  damping  term,  was  positive  where 
a  negative  value  was  expected.  It  was  found  that  Cn^  and  Cn^  were  linked;  if  one  was  fixed 
while  the  other  was  identified,  the  sum  of  the  two  effects  was  the  same  as  if  both  were  being 
identified.  Hence  the  ^  terms  were  fixed  at  zero  and  their  contribution  was  accounted  for  by 
identifying  the  yaw  rate  derivatives,  i.e.  Cy^,  Ci^  and  C„^.  The  match  was  still  good,  and  Cn^ 
had  the  correct  sign  for  all  cases  and  reasonable  magnitudes,  although  it  does  exhibit  a  large 
degree  of  scatter  about  the  wind  tunnel  data  points.  This  is  partly  a  result  of  the  aircraft  con¬ 
trol  system  being  engaged  during  the  flight  test  programme.  The  large  rapid  control  inputs, 
necessary  for  the  identification  process,  are  quickly  damped  by  the  control  system,  resulting  in 
reduced  identification  accuracy  of  the  yaw  damping  derivative, 

For  the  16°  wing  sweep  case  the  flight  test  results  agree  well  with  the  wind  tunnel  data.  These 
trends  are  reinforced  by  similar  trends  for  the  clean  configuration  in  [20],  [21]  and  [22].  Results 
for  the  roll  damping  derivative  Ci^,,  while  of  the  correct  sign,  are  approximately  25%  lower  in 
magnitude  than  the  wind  turmel  results  for  the  26°  wing  sweep  case,  but  display  similar  trends. 


5.2.3  Control  Derivatives 

The  lateral  control  derivatives  consist  of  side  force,  rolling  and  yawing  moments  due  to  aileron, 
rudder  and  spoiler  deflections,  Cy^^ ,  Cy^^ ,  Cy^^^ ,  and  ■ 

For  wing  sweeps  of  less  than  47°,  roll  control  activates  both  aileron  and  spoiler.  The  moments 
produced  by  these  controls  are  similar  and  so  their  individual  contributions  are  difficult  to  sep¬ 
arate.  Neither  spoilers  nor  ailerons  contribute  significantly  to  the  total  side  force,  and  ailerons 
cause  little  yawing  moment.  Accordingly,  the  derivatives  Cy^^ ,  Cy^^^  and  are  of  minor  im¬ 
portance.  In  order  to  improve  the  identification  process  they  were  therefore  fixed  at  their  a  priori 
values  and  only  Cy^^ ,  ,  Cig^^,Cns^  and  were  identified. 

The  non-dimensional  sideforce  due  to  rudder  deflection  derivative,  Cy^^  is  small  in  magnitude 
and  of  secondary  importance.  The  flight  test  results,  in  general,  agree  well  with  the  wind  hm- 
nel  data  for  16°  wing  sweep.  At  26°  wing  sweep,  the  flight  test  results  display  considerably 
more  scatter  about  the  wind  tunnel  data.  The  non-dimensional  rolling  moment  due  to  rudder 
deflection  derivative,  is  also  small  in  magnitude  and  of  secondary  importance.  The  flight 
test  results  are  scattered  about  the  wind  tunnel  data  at  16°  wing  sweep.  At  26°  wing  sweep  the 
flight  test  results  agree  well  with  the  wind  tunnel  data. 

The  non-dimensional  rolling  moment  due  to  aileron  deflection  derivative,  shows  consider¬ 
able  scatter  for  16°  wing  sweep  but  at  26°  wing  sweep  the  flight  test  results  agree  well  with  the 
wind  tunnel  data. 

The  non-dimensional  rolling  moment  due  to  spoiler  deflection  derivative,  is  a  primary  roll 
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control  derivative.  For  16°  wing  sweep  the  flight  test  results  show  good  agreement  with  the 
wind  tunnel  data.  At  26°  wing  sweep  the  results  are  consistently  lower  in  magnitude  than  the 
wind  tunnel  data  by  approximately  50%,  with  the  exception  of  a  single  point  at  an  angle-of- 
attack  of  9°. 

The  non-dimensional  yawing  moment  due  to  rudder  deflection  derivative,  Cng^  is  a  primary 
lateral  control  derivative.  For  both  wing  sweep  angles  the  flight  test  results  are  lower  in  magni¬ 
tude  than  the  wind  tunnel  data  by  as  much  as  50%.  This  is  similar  to  the  clean  aircraft  case,  [20], 
[21],  [22],  where  the  estimation  procedure  consistently  estimated  to  be  approximately  65% 

lower  in  magnitude  than  the  wind  tunnel  data  whilst  still  following  the  same  trends. 

The  non-dimensional  yawing  moment  due  to  spoiler  deflection  derivative,  is  a  small,  sec¬ 
ondary  lateral  control  derivative  and  the  flight  test  results  in  general  show  reasonable  agree¬ 
ment  with  the  wind  tunnel  data.  At  26°  wing  sweep  the  results  are  consistently  lower  in  mag¬ 
nitude  than  the  wind  tunnel  data  by  approximately  50%. 

It  should  be  noted  that  for  both  the  16°  and  26°  wing  sweep  cases  there  are  flight  test  points  that 
show  a  positive  ■  Positive  spoiler  deflection  corresponds  to  the  port  spoiler  deflected  up¬ 
wards,  as  is  shown  in  figure  1,  and  the  increased  drag  produced  by  the  port  wing  would  result 
in  a  negative  yawing  moment.  Therefore  the  derivative  is  expected  to  be  negative.  How¬ 
ever,  at  low  Mach  number  and  moderate  wing  sweep  the  derivative  is  reduced  and  any 
errors  introduced  into  the  estimation  of  become  significant,  the  noticeably  larger  Cramer- 
Rao  botmds  on  the  points  in  question  indicate  the  reduced  degree  of  confidence  in  these  results. 
The  reversed  sign  of  some  of  the  points  is  indicative  of  the  difficulty  involved  in  obtaining 

the  F-lllC  spoiler  derivatives.  Reference  [27]  details  additional  research  being  carried  out  into 
spoiler  characteristics. 


5.3  Summary 

The  flight-estimated  aerodynamic  derivatives  are  generally  in  agreement  with  the  wind  tunnel 
values.  There  are  three  exceptions;  the  pitch  stiffness,  Cma/  the  yaw  stiffness,  and  the  yaw¬ 
ing  moment  due  to  rudder  derivative,  .  The  flight  estimates  are  consistently  less  than  the 
wind  tunnel  derivatives  and  this  result  was  also  found  for  the  clean  aircraft  configuration,  [15], 
[16],  [17],  [20],  [21]  and  [22].  NASA  Dryden  reports  also  agree  with  this  trend  [26].  The  wind 
tunnel  results  for  the  yawing  moment  derivatives  of  the  F-lllC  are  not  unlike  those  of  other  air¬ 
craft  test  programs  carried  out  in  the  1960's.  These  derivatives  are  susceptible  to  model  sting- 
support  interference.  However,  no  details  of  the  wind  tunnel  facility,  model  size  and  accuracy 
or  the  test  procedure  are  available  to  confirm  this. 

The  augmented  pitch  and  yaw  damping  provided  by  the  automatic  flight  control  system,  whilst 
modifying  the  aircraft's  natural  response,  also  reduces  the  identifiability  of  a  number  of  sec¬ 
ondary  derivatives.  This  also  appears  to  have  contributed  to  the  reduction  in  the  values  of  the 
pitch  and  yaw  stiffness  derivatives  identified. 

The  results  also  indicate  that  for  the  26°  wing  sweep  case  there  is  a  tendency  to  underestimate 
the  spoiler  derivatives,  and  as  well  as  the  dihedral  effect,  Ci^.  The  clean  aircraft  con¬ 
figuration  does  not  show  this  trend.  There  is  no  obvious  reason  for  this  behaviour  and  the  NASA 
results  published  do  not  cover  the  take-off  and  landing  configurations. 
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Decoupling  the  spoiler  and  differential  stabilator,  i.e.  aileron,  derivatives  resulted  in  identifica¬ 
tion  difficulties.  Combining  the  spoiler  and  aileron  effects  [26],  whilst  improving  identifiability, 
does  not  allow  simple  comparison  with  wind  tunnel  data. 


6  CONCLUSIONS 

Stability  and  control  derivatives  describing  the  longitudinal  and  lateral  aerodynamic  character¬ 
istics  of  the  F-lllC  have  been  determined  from  flight  test  measurements.  The  results  presented 
in  this  report  are  for  16°  and  26°  wing  sweep  angles  with  the  aircraft  in  the  take-off  and  landing 
configurations.  Details  of  the  flight  test  schedule  and  the  data  processing  and  analysis  proce¬ 
dures  have  been  presented. 

The  flight  test  program  covered  only  limited  take-off  and  landing  test  points  thus  only  limited 
validation  over  the  wider  angle-of-attack  range  of  the  wind  tunnel  data  was  possible.  In  general 
all  primary  longitudinal  and  lateral  derivatives  show  satisfactory  repeatability  and  indicated 
similar  trends  to  the  wind  tunnel  data. 

Estimates  for  the  Cm^.Cnp  and  derivatives  were  consistently  lower  in  magnitude  than  the 
corresponding  wind  tunnel  data.  This  is  consistent  with  the  clean  aircraft  results  given  in  [15], 
[16],  [17],  [20],  [21]  and  [22].  For  the  26°  wing  sweep  cases,  the  two  spoiler  derivatives 
and  Cns^  ,  and  Ci^  were  underestimated.  Further  research  is  underway  into  obtaining  more 
accurate  estimates  of  spoiler  characteristics  [27]. 
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APPENDIX  A 

CALIBRATION  OF  ANGLE-OF-ATTACK  SENSORS 

The  angle-of-attack  sensor  position  errors  for  the  CADS  used  in  Phase  1,  and  NBTU  used  in 
Phase  3,  were  calculated  using  a  flight  path  reconstruction  computer  program  [28].  The  pro¬ 
gram  uses  measured  values  of  longitudinal  and  normal  acceleration  and  pitch  rate  as  inputs  to 
a  set  of  equations  describing  the  longitudinal  kinematic  relationships.  From  the  state  variables, 
the  body  axes  velocities  (v,  and  w)  and  pitch  attitude  (0),  the  output  quantities  angle-of-attack 
and  altitude  (a  and  H)  are  calculated.  The  flight  path  reconstruction  technique  uses  an  extended 
Kalman  filter  to  determine  the  calibration  parameters  relating  the  computed  and  measiued  out¬ 
put  variables.  For  the  angle-of-attack  instrument,  the  calibration  is  assumed  to  be  linear  with  a 
slope  Ka  and  bias  &«.  The  bias  term  is  not  needed  in  the  parameter  estimation  stage  since  bias 
errors  are  accounted  for  separately. 

The  flight  path  reconstruction  software  assumes  that  the  measured  angle  of  attack  may  be  de¬ 
scribed  by  an  equation  of  the  form : 


Oim  —  (f  T  •^a)®true  T  (■'^1) 

where  and  are  the  scale  factor  and  bias  respectively  and  the  subscripts  m  and  t  indicate  the 
measured  and  true  values  respectively.  In  the  flight  path  reconstruction  software,  terms  which 
adjust  the  a  measurement  to  the  eg  position  are  accoimted  for,  as  are  the  effects  of  the  vane  sen¬ 
sor  displacement  from  the  eg,  by  making  use  of  the  instrument  locations  with  respect  to  the  eg. 

The  aerodynamic  estimation  software,  pEst  [14],  defines  an  upwash  factor,  or  flow  amplifica¬ 
tion,  for  the  angle-of-attack  sensor,  Ka,  such  that : 

Q  P 

CX.Z  —  Koi\p-true\  ~  ”  ^cg)^  "F  (j/a  Veg)^  T  Olfe  (-^2) 

where  a^,  is  the  computed  response  or  the  equivalent  of  the  true  a. 

The  upwash  factor  in  pEst  can  be  expressed  in  a  similar  form  to  that  of  equation  A1 : 


=  +  (A3) 

The  variation  of  versus  true  airspeed  is  given  in  figure  A1  for  both  Phase  1  and  3  cases. 
While  there  is  some  variation  in  with  Mach  number  and  wing  sweep  angle,  these  variations 
are  not  significant  within  the  accuracy  of  the  data  and  so  a  constant  value  of  1.04  was  chosen  for 
Phase  1  compared  with  0.94  for  the  cruise  configuration.  An  investigation  into  the  effect  of 
variation  showed  that  errors  in  of  10  %  resulted  in  an  adjustment  of  typically  5  %  in  the  major 

derivatives.  The  difference  between  the  clean  and  take-off  and  landing  configuration  results  is 
postulated  to  be  a  consequence  of  the  upstream  influence  of  the  aerodynamically  'dirty'  take-off 
and  landing  configuration. 

Only  three  Phase  3  cases  were  available.  Hence  the  value  calculated  for  each  case  was  used. 
Koc  values  of  approximately  1.17  were  calculated  for  the  Phase  3  take-off  and  landing  configu¬ 
ration  cases  compared  with  a  value  of  1.06  for  the  clean  configuration. 
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APPENDIX  B 

CALIBRATION  OF  ANGLE-OF-SIDESLIP  SENSOR 


The  angle-of-sideslip  sensor  position  errors  for  the  CADS  used  in  Phase  1,  and  NBTU  used  in 
Phase  3,  were  calculated  using  the  full  six  degree-of-freedom  mode  of  the  flight  path  reconstruc¬ 
tion  computer  program  [28].  The  program  uses  measured  values  of  longitudinal,  lateral  and 
normal  accelerations  and  pitch,  roll  and  yaw  rates  as  inputs  to  a  set  of  equations  describing  the 
full  six  degree-of-freedom  kinematic  relationships.  From  the  state  variables,  the  body  axes  ve¬ 
locities  (-u,  V  and  w)  and  roll,  pitch  and  yaw  attitudes  6  and  V')/  the  output  quantities  of  ve¬ 
locity,  angle-of-attack,  angle-of-sideslip,  bank  and  pitch  angle  and  altitude  {V,  a,f3,4>,6  and  H) 
are  calculated.  The  flight  path  reconstruction  technique  uses  an  extended  Kalman  filter  to  deter¬ 
mine  the  calibration  parameters  relating  the  computed  and  measured  output  variables.  For  the 
angle-of-sideslip  instrument,  the  calibration  is  assumed  to  be  linear  with  a  slope  Kjs  and  bias  bp. 
The  bias  term  is  not  needed  in  the  parameter  estimation  stage  since  bias  errors  are  accounted  for 
separately. 

The  flight  path  reconstruction  software  assumes  that  the  measured  angle  of  sideslip  may  be  de¬ 
scribed  by  an  equation  of  the  form  : 


Pm  —  (1  +  ^fi)Ptrue  +  bp  (Bl) 

where  Xp  and  bp  are  the  scale  factor  and  bias  respectively  and  the  subscripts  m  and  t  indicate  the 
measured  and  true  values  respectively.  In  the  flight  path  reconstruction  software,  terms  which 
adjust  the  P  measurement  to  the  eg  position  are  accounted  for,  as  are  the  effects  of  the  vane  sen¬ 
sor  displacement  from  the  eg,  by  making  use  of  the  instrument  locations  with  respect  to  the  eg. 

The  aerodynamic  estimation  software,  pEst  [14],  defines  a  sidewash  factor,  or  flow  amplifica¬ 
tion,  for  the  angle-of-sideslip  sensor,  Kp,  such  that : 

Pz  =  Kp[Ptrue]  +  [zp  -  Zeg)^  -  {xp  -  X^g)^  +  Pb  (B2) 

where  P^  is  the  computed  response  or  the  equivalent  of  the  true  p. 

The  sidewash  factor  in  pEst  can  be  expressed  in  a  similar  form  to  that  of  equation  Bl : 


=  (1  +  A^)  (B3) 

The  variation  of  Kp  with  true  airspeed  is  given  in  figure  B2  for  both  Phase  1  and  Phase  3  cases. 
The  CADS  angle-of-sideslip  sensor  overestimated  the  true  value  by  approximately  55%  for  Phase 
1  take-off  and  landing  configuration  cases.  The  clean  aircraft  configuration  results  indicated  that 
Kp  was  overestimated  by  approximately  60%.  As  with  the  angle-of-attack  sensor  only  three 
cases  were  available  from  the  Phase  3  flights  hence  the  individual  Kp  values  were  used  in  the 
analysis.  The  NBTU,  as  expected,  gave  closer  estimates  of  the  true  angle-of-sideslip,  varying  by 
only  ±10%.  This  result  is  similar  to  that  obtained  for  the  clean  configuration  case. 

An  investigation  into  the  affect  of  Kp  variation  showed  that,  in  general,  changes  in  Kp  of  10  % 
resulted  in  maximum  changes  of  5  %  in  the  derivatives. 
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APPENDIX  C  -  SAMPLE  LONGITUDINAL  pEst  OUTPUT 

pEst:  Slanted  text  represents  user  inputs.  Units  are  defined  in  pEst  manual  [14]. 


pEst  program  for  parameter  estimation 
Richard  Maine/James  Murray  -  NASA  Dryden 
version  2.3.4  10  Aug  90 
this  run  date:  time: 

Help  is  available 


pEst:  read  m.plfSeOlc 

reading  measured  time  history  from  file  "m.plf6e01c" 
signal  dl  not  found.  Using  0. 
signal  d2  not  foimd.  Using  0. 
signal  d4  not  found.  Using  0. 
signal  thrust  not  found.  Using  0. 

pEst:  restore  c.plfSeOlc.final 

reading  program  status  from  file  "c.plf6e01c.finar' 

checking  consistency  of  file  "c.plf6e01c.finar'  and  file  "m.plfbeOlc" 

pEst:  states 

name  on?  limit  name  on?  limit  name  on?  limit 
alpha  T  .lOOOE+05  q  T  .lOOOE+05  theta  T  .lOOOE+05 
pEst:  resp 

name  weight  on?  name  weight  on? 

alpha  27.00  T  theta  13.00  T 
q  13.00  T  an  200.0  T 

pEst:  sho  stats 
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signal 

average 

signal 

average 

signal 

average 

de 

-1.58 

q 

-.145E-02 

ay 

-.198E-01 

dl 

.OOOE+00 

theta 

4.21 

pdot 

-.824 

d2 

.OOOE+00 

an 

1.00 

rdot 

-1.02 

da 

.845E-01 

ax 

-.883E-01 

qbar 

193. 

dr 

-.425 

qdot 

1.14 

mach 

.363 

d3 

-.224 

beta 

.339 

thrust 

.OOOE+00 

d4 

.OOOE+00 

P 

-.386 

alt 

353. 

V 

405. 

r 

-.650E-01 

alpha 

5.97 

phi 

-1.08 

pEst:  do  mmle3:para.plf6e01c 
reading  from  command  file. 

pEst:  par  xa  -25.4766 
name  value  on? 

xa  -25.48  F 

pEst:  par  ya  0.0000 
name  value  on? 

ya  .OOOOE+00  F 

pEst:  par  za  -1.5204 
name  value  on? 

za  -1.520  F 

pEst:  par  xb -32.1599 
name  value  on? 

xb  -32.16  F 

pEst:  paryb  0.0000 
name  value  on? 

yb  .OOOOE+00  F 

pEst:  par  zb  -3.5121 
name  value  on? 

zb  -3.512  F 

pEst:  par  xax  -2.5416 
name  value  on? 
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xax  -2.542  F 

pEst:  par  yax  -2.4308 
name  value  on? 

yax  -2.431  F 

pEst:  par  zax  -1.4371 
name  value  on? 

zax  -1.437  F 

pEst:  par  xay  -2.4391 
name  value  on? 

xay  -2.439  F 

pEst:  par  yay  -2.4583 
name  value  on? 

yay  -2.458  F 

pEst:  par  zay -1.4371 
name  value  on? 

zay  -1.437  F 

pEst:  parxan  -2.5416 
name  value  on? 

xan  -2.542  F 

pEst:  par  yan  -2.3617 
name  value  on? 

yan  -2.362  F 

pEst:  par  zan  -1.4371 
name  value  on? 

zan  -1.437  F 

pEst:  const  ix  73674.3359 
name  value 

ix  .7367E+05 
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pEst:  const  iy  387397.3750 
name  value 

iy  .3874E+06 

pEst:  const  iz  453923.0313 
name  value 

iz  .4539E+06 

pEst:  const  ixz  2542.6226 
name  value 

ixz  2543. 

pEst:  const  mass  2334.9495 
name  value 

mass  2335. 

pEst:  const  xcg  0.0 
name  value 

xcg  .OOOOE+00 

pEst:  const  ycg  0.0 
name  value 

ycg  .OOOOE+00 

pEst;  const  zcg  0.0 
name  value 

zcg  .OOOOE+00 

pEst:  resp  alpha  w=18 
name  weight  on? 

alpha  18.00  T 

pEst:  resp  q  w=13 
name  weight  on? 

q  13.00  T 

pEst:  resp  theta  w=13 
name  weight  on? 
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theta  13.00  T 

pEst:  resp  an  w=200 
name  weight  on? 

an  200.0  T 

pEst:  par  cnormO 
name  value  on? 

cNormO  .1767  T 

pEst:  par  cnorma 
name  value  on? 

cNorma  .1025  T 

pEst:  par  cnormde 
name  value  on? 

cNormde  .1543E-01  T 

pEst:  par  cma 
name  value  on? 

cma  -.8652E-02  T 

pEst;  par  cmq 
name  value  on? 

cmq  -35.36  T 

pEst; 

command  file  done. 
pEst:  it  20 


name 

value 

on? 

name 

value 

on? 

name 

value 

on? 

cNormO 

.1767 

T 

cNormde 

.1543E-01 

T 

cmq 

-35.36 

T 

cNorma 

.1025 

T 

cma 

-.8652E-02 

T 

cmde 

-.3059E-01 

T 

total  cost  =  59.15 
cost  per  response  : 
alpha  q  theta  an 
56.9  38.4  29.2  112. 

Iteration  1 

Lev-Marq  used  1  trials,  v  =  .OOOOE+00  len  =  .3121E-01  cost  =  59.15 
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name 

value 

on? 

delta 

name 

value 

on? 

delta 

cNormO 

.1767 

T 

.6425E-04 

cma 

-.8652E-02 

T 

-.3889E-06 

cNorma 

.1025 

T 

-.6514E-05 

cmq 

-35.36 

T 

.7772E-02 

cNormde 

.1543E-01 

T 

-.2783E-06 

cmde 

-.3059E-01 

T 

.1442E-06 

total  cost  =  59.15 


cost  per  response  : 


alpha  q  theta 
56.9  38.4  29.2 


an 

112. 


iteration  1  used  10  integrations 


**  iteration  converged 
estimate  used  10  integrations 
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pEst:  par  all 


name 

value 

on? 

name 

value 

on? 

name 

value 

cNormO 

.1767 

T 

cmd2 

•OOOOE+OO 

F 

qBias 

.OOOOE+OO 

cNorma 

.1025 

T 

cmd3 

.OOOOE+00 

F 

thetaBias 

.OOOOE+OO 

cNorma2 

.OOOOE+00 

F 

caO 

.OOOOE+OO 

F 

anBias 

.OOOOE+OO 

cNormAdot 

2.544 

F 

caa 

.OOOOE+OO 

F 

axBias 

.OOOOE+OO 

cNormq 

5.806 

F 

caa2 

.OOOOE+OO 

F 

qdotBias 

.OOOOE+OO 

cNormde 

.1543E-01 

T 

caAdot 

.OOOOE+OO 

F 

ka 

1.040 

cNormde2 

.OOOOE+00 

F 

caq 

.OOOOE+OO 

F 

xa 

-25.48 

cNormdl 

.OOOOE+00 

F 

cade 

.OOOOE+OO 

F 

ya 

.OOOOE+OO 

cNormd2 

.OOOOE+00 

F 

cade2 

.OOOOE+OO 

F 

za 

-1.520 

cNormdS 

.OOOOE+00 

F 

cadi 

.OOOOE+OO 

F 

xan 

-2.542 

cmO 

.OOOOE+00 

F 

cad2 

.OOOOE+OO 

F 

yan 

-2.362 

cma 

-.8652E-02 

T 

cad3 

.OOOOE+OO 

F 

zan 

-1.437 

cma2 

.OOOOE+00 

F 

vO 

.OOOOE+OO 

F 

xax 

-2.542 

cmAdot 

-4.665 

F 

alphaO 

.OOOOE+OO 

F 

yax 

-2.431 

cmq 

-35.36 

T 

qO 

.OOOOE+OO 

F 

zax 

-1.437 

cmde 

-.3059E-01 

T 

thetaO 

.OOOOE+OO 

F 

XV 

.OOOOE+OO 

cmde2 

.OOOOE+00 

F 

vBias 

.OOOOE+OO 

F 

yv 

.OOOOE+OO 

cmdl 

.OOOOE+00 

F 

alphaBias 

.OOOOE+OO 

F 

zv 

.OOOOE+OO 

name 

value 

on? 

name 

value 

on? 

name 

value 

on? 

cyO 

.7402E-01 

F 

cld3 

-.1600E-02 

F 

rBias 

.OOOOE+OO 

F 

cyb 

-.1230E-01 

F 

cld4 

.OOOOE+OO 

F 

phiBias 

.OOOOE+OO 

F 

cyb3 

.OOOOE+OO 

F 

cnO 

.7418E-03 

F 

ayBias 

-.5902E-01 

F 

cyBdot 

.OOOOE+OO 

F 

cnb 

.2135E-02 

F 

pdotBias 

.OOOOE+OO 

F 

cyp 

.OOOOE+OO 

F 

cnb3 

.OOOOE+OO 

F 

rdotBias 

.OOOOE+OO 

F 

cyr 

.OOOOE+OO 

F 

cnBdot 

.OOOOE+OO 

F 

kb 

1.640 

F 

cyda 

.7661E-03 

F 

cnp 

-.5000E-01 

F 

xb 

-32.16 

F 

cydr 

.3829E-02 

F 

cnr 

-.2172 

F 

yb 

.OOOOE+OO 

F 

cyd3 

-.8200E-04 

F 

end  a 

.2827E-02 

F 

zb 

-3.512 

F 

cyd4 

.OOOOE+OO 

F 

endr 

-.3189E-02 

F 

xay 

-2.439 

F 

clO 

.8894E-02 

F 

cnd3 

-.3800E-03 

F 

yay 

-2.458 

F 

clb 

-.1403E-02 

F 

cnd4 

.OOOOE+OO 

F 

zay 

-1.437 

F 

clb3 

.OOOOE+OO 

F 

betaO 

.OOOOE+OO 

F 

gAlpha 

.OOOOE+OO 

F 

clBdot 

.OOOOE+OO 

F 

pO 

.OOOOE+OO 

F 

gQ 

.OOOOE+OO 

F 

clp 

-.4295 

F 

rO 

.OOOOE+OO 

F 

gBeta 

.OOOOE+OO 

F 

clr 

.1000 

F 

phiO 

.OOOOE+OO 

F 

gP 

.OOOOE+OO 

F 

clda 

.2405E-02 

F 

betaBias 

.OOOOE+OO 

F 

gP 

.OOOOE+OO 

F 

cldr 

-.1005E-02 

F 

pBias 

.OOOOE+OO 

F 

on? 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 


F 

F 
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APPENDIX  D  -  SAMPLE  LATERAL  pEst  OUTPUT 

pEst:Sianted  fexi  represents  user  inputs.  Units  are  defined  in  pEst  manual  [14]. 


pEst  program  for  parameter  estimation 
Richard  Maine /James  Murray  -  NASA  Dryden 
version  2.3.4  10  Aug  90 
this  run  date:  time: 

Help  is  available 


pEst:  read  m.plf6e38c 

reading  measured  time  history  from  file  "m.plf6e38c" 
signal  dl  not  found.  Using  0. 
signal  d2  not  foimd.  Using  0. 
signal  d4  not  found.  Using  0. 
signal  thrust  not  found.  Using  0. 

pEst:  restore  c.plf6e38c.£inal 

reading  program  status  from  file  "c.plf6e38c.fmal" 

checking  consistency  of  file  "c.plf6e38c.final"  and  file  "m.plf6e38c" 


pEst:  states 
name  on? 

limit 

name 

on? 

limit 

beta  T 

.lOOOE+05 

r 

T 

.lOOOE+05 

P  T 

.lOOOE+05 

phi 

T 

.lOOOE+05 

pEst:  resp 

name 

weight 

on? 

name 

weight 

on?  name 

weight 

on? 

beta 

7.500 

T 

r 

40.00 

T  ay 

1333. 

T 

P 

4.400 

T 

phi 

4.000 

T 

pEst:  sho  stats 
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signal 

average 

signal 

average 

signal 

average 

de 

-5.62 

q 

.223 

ay 

-.146E-01 

dl 

.OOOE+00 

theta 

8.14 

pdot 

-.268 

d2 

.OOOE+00 

an 

.990 

rdot 

-.807 

da 

-.452E-01 

ax 

-.154 

qbar 

164. 

dr 

-.225 

qdot 

1.16 

mach 

.340 

d3 

.332 

beta 

.662 

thrust 

.OOOE+00 

d4 

.OOOE+00 

P 

-.597 

alt 

.115E+04 

V 

378. 

r 

-.172 

alpha 

10.4 

phi 

-.671 

pEst:  do  mmleS.para 
reading  from  command  file. 

pEst:  par  xa  -25.5547 
name  value  on? 

xa  -25.55  F 

pEst:  parya  .0000 
name  value  on? 

ya  .OOOOE+00  F 

pEst:  par  za  -1.5037 
name  value  on? 

za  -1.504F 

pEst:  par  xb  -32.2380 
name  value  on? 

xb  -32.24  F 

pEst:  paryb  .0000 
name  value  on? 

yb  .OOOOE+00  F 

pEst:  par  zb  -3.4954 
name  value  on? 

zb  -3.495  F 

pEst:  par  xax  -2.6197 
name  value  on? 
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xax  -2.620  F 

pEst:  par  yax  -2.4308 
name  value  on? 

yax  -2.431  F 

pEst:  par  zax  -1.4204 
name  value  on? 

zax  -1.420  F 

pEst:  par  xay  -2.5172 
name  value  on? 

xay  -2.517  F 

pEst:  par  yay  -2.4583 
name  value  on? 

yay  -2.458  F 

pEst:  par  zay  -1.4204 
name  value  on? 

zay  -1.420  F 

pEst:  par  xan  -2.6197 
name  value  on? 

xan  -2.620  F 

pEst:  paryan  -2.3617 
name  value  on? 

yan  -2.362  F 

pEst:  par  zan  -1.4204 
name  value  on? 

zan  -1.420  F 

pEst:  const  ix  69022.4453 
name  value 

ix  .6902E+05 
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pEst:  const  iy  365324.0625 
name  value 

iy  .3653E+06 

pEst:  const  iz  427141.6875 
name  value 

iz  .4271E+06 

pEst:  const  ixz  4071.1189 
name  value 

ixz  4071. 

pEst:  const  mass  2250.4507 
name  value 

mass  2250. 

pEst:  const  xcg  0.0 
name  value 

xcg  .OOOOE+00 

pEst:  const  ycg  0.0 
name  value 

ycg  .OOOOE+00 

pEst:  const  zcg  0.0 
name  value 

zcg  .OOOOE+00 

pEst: 

command  file  done. 
pEst:  it  20 


name 

value 

on? 

name 

value 

on? 

name 

value 

0] 

cyO 

.6118E-02 

T 

clp 

-.3053 

T 

cnb 

.9239E-03 

T 

cyb 

-.1737E-01 

T 

clda 

-.1239E-02 

T 

cnr 

-.1350 

T 

cydr 

.6585E-02 

T 

cldr 

.4348E-03 

T 

cndr 

-.1094E-02 

T 

clO 

.1007E-02 

T 

cld3 

-.4697E-03 

T 

cnd3 

-.7729E-04 

T 

clb 

-.1715E-02 

T 

cnO 

-.6007E-03 

T 

ayBias 

-.1430E-01 

T 

total  cost  =  45.80 
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cost  per  response : 
beta  p  r  phi  ay 

11.8  21.7  94.0  39.6  62.0 

Iteration  1 

Lev-Marq  used  1  trials,  v  =  .OOOOE+00  len  =  .1673  cost  =  45.80 


name 

value 

on? 

delta 

name 

value 

on? 

delta 

cyO 

.6136E-02 

T 

.1819E-04 

cld3 

-.4696E-03 

T 

.6128E-07 

cyb 

-.1737E-01 

T 

.1348E-05 

cnO 

-.6001E-03 

T 

.5786E-06 

cydr 

.6586E-02 

T 

.6603E-06 

cnb 

.9233E-03 

T 

-.5162E-06 

clO 

.1007E-02 

T 

.9192E-07 

cnr 

-.1347 

T 

.2231E-03 

clb 

-.1716E-02 

T 

-.7743E-06 

cndr 

-.1094E-02 

T 

.4461E-06 

clp 

-.3052 

T 

.9662E-04 

cnd3 

-.7735E-04 

T 

-.5351E-07 

clda 

-.1239E-02 

T 

-.4190E-06 

ayBias 

-.1433E-01 

T 

-.2753E-04 

cldr 

.4357E-03 

T 

.8805E-06 

total  cost  =  45.80 
cost  per  response : 

?.i  rr  0X1  iteration  1  used  1 9  integrations 

11.8  21.7  94.1  39.5  61.9  ^ 

**  iteration  converged 

estimate  used  19  integrations 
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pEst:  par  all 


name 

value 

on? 

name 

value 

on? 

name 

value 

cNormO 

.OOOOE+00 

F 

cmd2 

.OOOOE+OO 

F 

qBias 

.OOOOE+OO 

cNorma 

.1000 

F 

cmd3 

.OOOOE+OO 

F 

thetaBias 

.OOOOE+OO 

cNorma2 

.OOOOE+00 

F 

caO 

.OOOOE+OO 

F 

anBias 

.OOOOE+OO 

cNormAdot 

.OOOOE+00 

F 

caa 

.OOOOE+OO 

F 

axBias 

.OOOOE+OO 

cNormq 

•OOOOE+OO 

F 

caa2 

.OOOOE+OO 

F 

qdotBias 

.OOOOE+OO 

cNormde 

.5000E-02 

F 

caAdot 

.OOOOE+OO 

F 

ka 

1.000 

cNormde2 

.OOOOE+OO 

F 

caq 

.OOOOE+OO 

F 

xa 

-25.55 

cNormdl 

.OOOOE+OO 

F 

cade 

.OOOOE+OO 

F 

ya 

.OOOOE+OO 

cNormd2 

•OOOOE+OO 

F 

cade2 

.OOOOE+OO 

F 

za 

-1.504 

cNormdS  , 

.OOOOE+OO 

F 

cadi 

.OOOOE+OO 

F 

xan 

-2.620 

cmO 

.OOOOE+OO 

F 

cad2 

.OOOOE+OO 

F 

yan 

-2.362 

cma 

-.lOOOE-Ol 

F 

cad3 

.OOOOE+OO 

F 

zan 

-1.420 

cma2 

.OOOOE+OO 

F 

vO 

.OOOOE+OO 

F 

xax 

-2.620 

cmAdot 

.OOOOE+OO 

F 

alphaO 

.OOOOE+OO 

F 

yax 

-2.431 

cmq 

-20.00 

F 

qO 

.OOOOE+OO 

F 

zax 

-1.420 

cmde 

-.2000E-01 

F 

thetaO 

.OOOOE+OO 

F 

XV 

.OOOOE+OO 

cmde2 

.OOOOE+OO 

F 

vBias 

.OOOOE+OO 

F 

yv 

.OOOOE+OO 

cmdl 

.OOOOE+OO 

F 

alphaBias 

.OOOOE+OO 

F 

zv 

.OOOOE+OO 

name 

value 

on? 

name 

value 

on? 

name 

value 

on? 

cyO 

.6136E-02 

T 

cld3 

-.4696E-03 

T 

rBias 

.OOOOE+OO 

F 

cyb 

-.1737E-01 

T 

cld4 

.OOOOE+OO 

F 

phiBias 

.OOOOE+OO 

F 

cyb3 

•OOOOE+OO 

F 

cnO 

-.6001E-03 

T 

ayBias 

-.1433E-01 

T 

cyBdot 

•OOOOE+OO 

F 

cnb 

.9233E-03 

T 

pdotBias 

.OOOOE+OO 

F 

cyp 

.1884 

F 

cnb3 

.OOOOE+OO 

F 

rdotBias 

.OOOOE+OO 

F 

cyr 

.3013 

F 

cnBdot 

.OOOOE+OO 

F 

kb 

1.550 

F 

cyda 

-.7500E-03 

F 

cnp 

-.1956E-01 

F 

xb 

-32.24 

F 

cydr 

.6586E-02 

T 

cnr 

-.1347 

T 

yb 

.OOOOE+OO 

F 

cyd3 

-.8200E-04 

F 

cnda 

-.1312E-03 

F 

zb 

-3.495 

F 

cyd4 

.OOOOE+OO 

F 

cndr 

-.1094E-02 

T 

xay 

-2.517 

F 

clO 

.1007E-02 

T 

cnd3 

-.7735E-04 

T 

yay 

-2.458 

F 

clb 

-.1716E-02 

T 

cnd4 

.OOOOE+OO 

F 

zay 

-1.420 

F 

clb3 

.OOOOE+OO 

F 

betaO 

.OOOOE+OO 

F 

gAlpha 

.OOOOE+OO 

F 

clBdot 

.OOOOE+OO 

F 

pO 

.OOOOE+OO 

F 

gQ 

.OOOOE+OO 

F 

clp 

-.3052 

T 

rO 

.OOOOE+OO 

F 

gBeta 

.OOOOE+OO 

F 

clr 

.1499 

F 

phiO 

.OOOOE+OO 

F 

gP 

.OOOOE+OO 

F 

clda 

-.1239E-02 

T 

betaBias 

.OOOOE+OO 

F 

gR 

.OOOOE+OO 

F 

cldr 

.4357E-03 

T 

pBias 

.OOOOE+OO 

F 

on? 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 


F 

F 
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APPENDIX  E  -  LONGITUDINAL  EQUATIONS  OF  MOTION 


State  vector 

X  =  (a.g,  6>) 

Control  vector 

u  =  (^e) 

Observation  vector 

^  —  (<^m ,  qm ,  ,  9m) 

The  nonlinear  longitudinal  state  equations  are  : 

qSR  gR 

q - +  —(cos  9  cos  a  +  sin  0  sin  a) 

mV  V  ’ 


a 


TR 


mV 

lyyq  =  qScRC„ 


■  sin  a 


=  ^0  +  q 


The  ^0  term  is  included  to  allow  for  instrument  biases. 


The  longitudinal  observation  equations  are  ; 


—  CKq  +  Ka{ 

Q 

1 

Qm 

=  qo  +  q 

Om 

=  6q-\- 6 

^rim 

,  qS 

—  Ono  H - 

CN+~q  + 

mg 

gR 

R^g 

-  n 

~  ^Xq 

mg 

CA  +  ^q- 
gR 

T 

mg 

qm 

=  90 +  g 

The  ao,  90,  ^0,  ,  ax^  and  go  terms  represent  instrument  biases. 


DSTO-TR-0321 


The  expansions  of  the  longitudinal  force  and  moment  coefficients  are 
Cn  =  Cno  +  +  Cn,  +  Cns^ 

Cm  —  Cmo  +  Cmc,Oi  +  Cm^-^^  + Cms^Se  +  Cmi, 
Ca  =  Cao  +CAc<^  +  CA^-^^  +  CAsJe 
Cl  =  Cn  cos  a  —  Ca  sin  a 


ac 

2VR 
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APPENDIX  F  -  LATERAL  EQUATIONS  OF  MOTION 

State  vector  x  =  (a,  (5,  p,  q,  r,  0,  <f>) 

Control  vector  u  =  (^a,  ^e,  ^sp) 

Observation  vector  Z  —  ( 1  Pm  t  Pm  1  Qm:  ‘^m  J  J  4^m )  ^Xm  T  ^Vtti  j  ) 

The  nonlinear  lateral-directional  state  equations  are  : 
a 

P 

pix  ~  ^ ^xz 

qly  ~  '^lyz  V^xy 
'i'lz  ~  pixz  ^lyz 

e 

The  $0  arid  c^o  are  included  to  allow  for  instrument  biases. 

The  lateral  observation  equations  are  : 


=  ao  + 

,  2/a  \ 

Pm 

=  Po  +  Kp{p- 

‘7*’  + V* 

Pm 

=  Po  +P 

9m 

=  90  +  9 

qSR  ^  gR  .  ^  ^  , 

1  ~  'hr _ a  „  (cOS  g  COS  0  COS  Q!  +  Sin  0  Sill  a) 


mV  cos  P  V  cos  P 
—  ta.n  P{pcosa  -f  rsin  a) 


TR 


mV  cos  P 


sin  a 


— i2(cos/3cos0sin  (f)  —  sin  /3 (cos  0  cos  iji)  sin  a  —  sin  ^cosa)) 
qS 

-\-p  sin  a  —  r  cos  a  -| - —RCy 

mV 


qShRCi  +  {qr(Iy  -  h)  +  (q^  -  r^)Iy^  +  pql^^  -  prI^y)/R 


-  qSbRCm  +  (prih  -  h)  +  (r  -  p  )Ixz  +  qrixy  -  pqIyz)/R 


qSbRCn  +  {pq{Ix  -  ly)  +  [p  -  q  )Ixy  +  priyz  -  qrla:z)lR 


=  6q  +  q  cos  0  —  r  sin 


=  ^0  +  P  +  tan  0(r  cos  -f  g  sin 
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Tm  -  ro  +  r 
Gm  =  ^0  +  ^ 
<Pm  —  ^0  +  ^ 


^X-m. 

—  ^Xo  ■“ 

^Oa  +  (%  - 
mg  gR 

Vax  •\ 
——r)  — 

gR  ’ 

^ax 

R^g 

(g2  +  r2)+  — 
mg 

O'Vm 

—  fflyo  + 

^Cy  +  - 

mg  gR 

^ay  ... 

Hay 

R^g 

(p2  +  r2) 

^rim 

—  O-no  T 

^Cm  +  9  - 

mg  gR 

.  + 

r_2  1  „2\ 

+^> 

Pm 

=  Po+P 

Qm 

=  qo  +  q 

Sm 

=  ro  +  r 

The  ao,Po,po,qo,ro,eo,<f>o,  dxo )  ^yo  > 


a„o  ,Vo,Qo  and  tq  terms  represent  instrument  biases. 


The  expansions  of  the  lateral  force  and  moment  coefficients  are  : 


pb 


rb 


Cy  =  Cyo  +  CYp/3  +  Cyp  +  CYr  +  Cy^S 


pb 


rb 


Cl  =  C/o  +  Ci„(3  -\-Cip +  Cl,  +  C//  +  Cl- 


13b 


Cn  —  Cno  +  CnnP -\- Cnp  Cn,  +  Cns^  +  Cn^ 


2YPI  --  2VR 


2VR 


where  the  5  terms  are  summed  over  all  controls. 


Cys  5  =  Cys^  +  Cys,  Sr  +  Cys^^  ^sp 

Cl^S  =  Cl^Ja  +  Cl^^Sr  +  Cl^^^Ssp 

Cng  5  =  Cng^  Sa  +  C„g,  Sr  +  Cfig^^  Sgp 
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